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1)  INTRODUCTION: 


We  have  developed  a  new  technique  for  obtaining  4D  dynamic  contrast-enhanced  (DCE)  digital 
breast  tomosynthesis  (DBT)  images  that  will  allow  the  measurement  of  tumor  perfusion  with 
improved  accuracy.  Current  contrast-enhanced  DBT  systems  produce  images  at  2-3  fixed  time 
points.  Our  novel  acquisition  technique  is  capable  of  producing  perfusion  images  at  30-60  time 
points  for  the  same  radiation  dose.  The  method  is  compatible  with  both  temporal  and  dual¬ 
energy  (DE)  subtraction  methods.  In  our  new  method,  a  single  projection  image  is  acquired,  for 
example,  once  every  5  seconds.  One  complete  tomosynthesis  projection  series  consists  of  a  set 
of  projection  images  acquired  at  distinct  angles.  In  the  new  method,  several  complete 
tomosynthesis  projection  series  are  acquired  sequentially.  Reconstruction  is  performed  using 
one  full  set  of  sequentially-acquired  images  (or  a  sufficient  subset);  however,  unlike 
conventional  DBT,  the  starting  angle  (and  hence  measurement  time  point)  is  arbitrary  in  the  new 
method.  The  resulting  4D  data  set  consists  of  many  time -resolved  3D  functional  measurements 
of  tumor  perfusion,  offering  the  potential  for  superior  lesion  characterization  and  hence 
diagnostic  accuracy.  In  this  grant,  we  have  built  and  tested  the  system,  including  extensive 
characterizations  of  spatial  resolution  and  temporal  resolution. 

2)  KEYWORDS:  Contrast-enhanced  imaging,  x-ray  imaging,  digital  breast  tomosynthesis 
(DBT),  4D  imaging,  modulation  transfer  function  (MTF),  spatial  resolution,  temporal 
resolution,  contrast  agent  quantification. 
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3)  OVERALL  PROJECT  SUMMARY: 


In  this  final  report,  we  can  demonstrate  excellent  progress  on  this  project.  We  presented  some 
aspects  of  this  research1  at  the  RSNA  Annual  Meeting  in  Chicago  in  November  2011,  and  will 
present  our  completed  results  at  the  RSNA  Annual  Meeting  in  Chicago  in  November  2013.  In 
addition,  we  have  one  manuscript  submitted  and  another  in  preparation. 

The  primary  purpose  of  this  grant  is  to  demonstrate  that  dynamic  4D  DCE-DBT  is 
technologically  feasible  and  results  in  superior  measurements  of  tumor  perfusion.  Two 
objectives  were  proposed:  (#1)  Build  a  DBT  prototype  to  demonstrate,  in  proof-of-principle, 
that  4D  DCE-DBT  is  technically  feasible.  (#2)  Develop  techniques  to  derive  accurate  measures 
of  breast  tumor  perfusion  dynamics.  We  have  successfully  completed  both  aims  and  are 
currently  in  the  process  of  preparing  a  follow-up  grant  application  to  bring  this  research  to  the 
next  logical  stage  -  the  development  of  a  clinical  prototype.  We  are  also  in  the  process  of 
preparing  a  patent  disclosure. 

Aim  1 :  System  Design  and  Fabrication 

A  conceptual  diagram  of  the  imaging  method  is  shown  in  Figure  1.  Rather  than  making  three 
rapid  acquisitions  of  13  projections  each,  thereby  giving  three  discrete  time  points  for  which  to 
perform  contrast-enhanced  imaging,  we  acquire  the  same  number  of  images  (and  thus  use  the 
same  dose);  however,  each  projection  image  is  delayed  from  the  next  by  a  finite  period  of  time. 
For  example,  in  the  example  shown,  a  delay  of  5  seconds  between  each  projection  image  would 
result  in  a  total  acquisition  time  of  190  sec  [5  sec  *  (3*13  -  1)].  Reconstructions  would  then  be 
made  at  27  discrete  time  points  [2*13  +  1].  In  this  aim,  we  proposed  to  build  a  proof-of- 
principle  prototype  for  both  the  acquisition  and  reconstruction  systems. 


Reconstructions 


Acquisitions 


Time 


Figure  1:  Illustration  of  the  acquisition  and  reconstruction  method 


We  have  completed  both  the  design  and  the  fabrication  of  the  imaging  system.  Projection 
images  are  acquired  using  a  custom  x  ray  imaging  testing  station  with  x  ray  source,  detector, 
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rotary  stage,  and  filter  as  shown  in  Figure  2.  The  testing  station  uses  a  GE  Senographe  500  x  ray 
generator  (General  Electric  Health  Care,  Chalfont-St.  Giles,  UK)  with  a  selenium  flat-panel 
detector  (Direct  Radiography  Corp,  Newark,  DE).  Parker  Motion  stepper  motors  (Parker 
Hannifin,  Rohnert  Park,  CA)  are  used  to  drive  a  rotary  stage  on  which  a  phantom  can  be 
mounted  and  rotated,  as  well  as  a  custom  filter  wheel;  magnetic  sensors  are  used  to  precisely 
control  the  position  of  the  rotary  stage  and  filter  wheel.  The  filter  wheel  is  laser-cut  acrylic  with 
various  mounted  filter  metals,  including  0.5  mm  of  aluminum  used  in  this  experiment. 
Additionally,  the  x  ray  tube  has  a  25  pm  molybdenum  filter  mounted  to  the  exit  port. 

In  a  real  clinical  system,  the  detector  and  x-ray  tube  would  move  while  the  patient  is  held 
stationary.  The  system  we  built  rotates  the  test  object,  this  approach  was  chosen,  as  it  is  easier  to 
make  while  remaining  mathematically  equivalent  to  the  clinical  approach. 


We  have  also  developed  a  method  of  preprocessing  an  existing  reconstruction  system  (Piccolo, 
Real-time  Tomography,  LLC)  available  in  our  lab  to  enable  reconstruction  of  the  resultant 
images. 


38.5  cm  87.5  cm 


test  tool 


Flat-panel 

detector 


Detector  Rotarystage  Collimatorand  X-raytube 

and  test  tool  filter  wheel 

Figure  2:  (top)  The  upper  image  depicts  a  schematic  of  custom  x-ray  imaging  station  from  aerial  view.  An  x-ray 
source,  aluminum  filter,  and  rotary  stage  with  test  object  are  positioned  in  succession  in  front  of  a  flat  panel  detector, 
(bottom)  The  image  below  is  a  photograph  of  the  x-ray  testing  station  (filter  wheel  is  not  visible). 


7 


Aim  2:  Quantitation  of  Iodine  Uptake 


We  have  developed  a  hybrid  image  subtraction  method  in  which  we  perform  temporal 
subtraction  of  the  dual-energy  images.  This  technique  has  several  advantages,  including  the 
ability  to  suppress  scatter  and  other  artifacts  in  the  dual  energy  images,  insensitivity  to  motion 
artifacts,  the  ability  to  quantify  contrast  agent  uptake,  and  good  noise  performance. 

The  ability  to  quantify  the  iodine  in  terms  of  areal  density  is  particularly  important.  The  data, 
shown  in  Figure  3,  illustrate  this  principle.  Shown  in  Figure  4a  is  the  subtracted  signal  as  a 
function  of  the  areal  density  of  the  contrast  material  in  a  static  phantom  already  in  use  in  the  lab2. 
The  measurements  are  breast  thickness  dependent.  Shown  are  the  results  for  breasts  2  cm  and 
9  cm  thick.  The  data  for  multiple  thicknesses  have  been  summarized  in  Fig  4b.  These  data 
allowed  us  to  quantify  the  uptake  of  tumors  in  absolute  numbers  for  the  first  time.  These  results 
were  presented  at  the  2011  Annual  Meeting  of  the  RSNA1. 


a)  b) 


Figure  3.  (a)  A  schematic  diagram  showing  one  layer  of  the  contrast  phantom.  The  inclusions  are  specified  in  Table 
1 .  (b)  A  hybrid  temporal  dual-energy  subtraction  radiograph  of  the  phantom  showing  the  iodine  inclusions. 
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Table  1:  The  phantom  illustrated  in  Figure  2  has  10 
iodinated  inclusions  with  the  specified  size  and  areal 
density. 


Lesion  depth 
[mm] 

Lesion  width 
[mm] 

Areal  density  [mg  I/cm2] 

5.035  mg  I/ml  2.575  mg  I/ml 

9.6 

10.6 

4.83 

2.47 

8 

8.8 

4.03 

2.06 

6.2 

7.0 

3.12 

1.60 

4.8 

5.4 

2.42 

1.24 

4 

4.6 

2.01 

1.03 

a) 


b) 


Areal  Density  of  Iodine  (mg/cm2) 


9 

8 

7 


£  6 
U 


1 


0 

0 


2  4  6  8  10 

Breast  Thickness  (cm) 


Figure  4.  (a)  Using  static  phantoms  of  different  thickness,  it  is  possible  to  quantify  the  signal  as  a  function  of  areal 
density.  Two  thicknesses  are  shown  for  illustration,  (b)  The  calibration  coefficient  (slope  of  lines  in  Fig  4a)  can 

then  be  calculated  for  each  breast  thickness. 
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It  was  necessary  to  estimate  the  uptake  accurately  as  a  function  of  time,  allowing  for  the  fact  that 
the  uptake  is  time  varying.  Through  experimentation,  we  have  found  that  the  uptake  is  given  by 
the  temporal  average  of  the  uptake  over  the  time  interval  used  in  the  reconstruction  of  the  image. 
In  the  current  reconstruction  scheme,  the  individual  source  projection  images  are  equally 
weighted  in  the  reconstruction  and  the  reconstruction  (over  the  key  steps)  is  linear.  As  such,  a 
moving  average  is  sufficient  to  model  the  temporal  properties  of  the  imaging  system. 

These  properties  were  tested  with  the  phantom  shown  in  Figure  5.  The  concentration  of  iodine 
in  the  circulating  fluid  was  controlled  by  the  addition  of  contrast  agent  or  distilled  water  to  the 
circulating  fluid.  Calibration  standards  with  various  concentrations  of  iodine  are  included  in  the 
phantom  to  act  as  a  test  standard. 


Figure  5:  A  flow  phantom  was  built  to  simulate  temporal  dependencies  of  blood  flow  in  the  breast.  The  phantom 
consists  of  a  main  flow  chamber  connected  to  a  pump  and  liquid  reservoir.  Contrast  agent  can  be  added  to  the 
reservoir  via  a  power  injector  (not  shown).  Six  calibration  chambers  are  aligned  along  the  top  edge  of  the  phantom. 


Figure  6:  A  sample  temporal  measurement  of  blood  flow  conducted  with  the  imaging  system.  The  actual  iodine 
concentration  is  shown  (black),  together  with  the  expected  value  (blue)  calculated  as  the  moving  average  of  the 
actual  data,  and  the  measured  results  (red).  Close  agreement  between  expected  and  observed  results  is  seen. 
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Figure  6  shows  the  results  of  one  such  temporal  concentration  measurement.  The  black  curve  is 
the  actual  concentration  as  detennined  by  the  addition  of  contrast  agent  and/or  distilled  water  to 
the  circulating  fluid.  The  blue  curve  is  the  expected  concentration  assuming  a  moving  average 
over  the  17  images  used  in  the  reconstruction.  The  red  curve  is  the  measured  concentration  as  a 
function  of  time.  It  can  be  seen  from  comparison  that  the  theoretical  and  experimental  data 
match  closely. 

During  the  above  experiments,  it  became  clear  that  it  was  possible  to  reconstruct  tomographic 
data  from  a  subset  of  angular  measurements;  for  example,  in  the  examples  above,  the  full  data 
acquisition  consisted  of  reconstructions  from  projection  images  measured  at  17  unique  angles. 
However,  it  is  possible  to  perfonn  the  reconstruction  with  any  number  of  unique  angles  (e.g.  9 
angles).  In  so  doing,  one  is  reducing  the  range  of  angles  used  (which  reduces  z-resolution),  and 
the  center  angle  will  vary  with  time. 

The  most  significant  effort  was  with  regard  to  the  spatial  resolution.  Given  that  this  problem  has 
similar  origins  to  the  problem  of  oblique  reconstructions  in  tomosynthesis3’ 4  (for  which  we  have 
recently  had  two  papers  accepted  to  Medical  Physics),  we  spent  considerable  time  and  effort  to 
understand  the  implications  of  these  changes. 

The  details  of  these  results  are  provided  in  the  attached  manuscript.  Provided  here  is  a  brief 
summary  of  the  results  seen  in  Figure  7.  The  MTF  is  computed  using  the  algorithm  described  by 
Carton  et  al5'  6 .  For  the  tomosynthesis  reconstructions,  the  MTF  is  computed  by  taking  the  MTF 
of  an  edge  in  a  slice  of  the  reconstruction.  Briefly,  the  procedure  is  as  follows: 

1 .  A  rectangular  region  of  interest  (of  750  x  350  pixels  or  80.3  x  37.5  mm)  including  the 
edge  is  selected. 

2.  A  derivative  filter  is  used  to  find  the  edge. 

3.  The  edge  is  interpolated  and  a  super-sampled  edge  spread  function  (ESF)  is  computed 
over  the  ROI  using  a  sliding  window  with  a  sub-region  column  size  of  20. 

4.  A  monotonic  filter  is  used  to  reduce  the  noise  of  the  ESF7. 

5.  The  line  spread  function  (LSF)  is  computed  by  differentiation  of  the  ESF  with  a 
differentiation  kernel  of  [-1  1]  and  correction  for  finite  element  differentiation. 

6.  The  MTF  is  computed  by  taking  the  Fourier  transform  of  the  LSF. 

The  effect  of  the  reconstruction  slice  offset  relative  to  the  plane  that  contains  the  edge  was 
measured  by  obtaining  reconstruction  slices  in  0. 1  mm  increments  from  the  true  slice  depth 
containing  the  edge.  This  was  perfonned  to  estimate  the  slice  spacing  that  is  necessary  to 
produce  a  reconstruction  data  set  that  is  minimally  degraded  with  respect  to  the  spatial 
resolution.  These  data  are  summarized  in  Figure  7B  and  1C.  It  was  found  that  the  optimal  slice 
spacing  is  dependent  upon  the  acquisition  geometry. 

The  effect  of  the  reconstruction  slice  angle  in  an  oblique  reconstruction  on  the  MTF  was 
assessed  by  changing  the  angle  of  the  lead  sheet  with  respect  to  the  detector  in  the 
reconstruction.  These  data  are  summarized  in  Figure  7D.  They  show  that  there  is  only  a  small 
reduction  in  spatial  resolution  as  the  angle  of  the  reconstructed  plane  is  varied  (e.g.,  in  oblique 
reconstructions). 
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Figure  7:  Various  characterizations  of  the  spatial  resolution  of  the  imaging  system.  A)  the  dependency  of  the  MTF 
on  the  acquisition  angular  range;  B)  the  dependency  of  the  resolution  on  the  plane  of  focus  of  the  reconstruction 
relative  to  the  position  of  the  edge  phantom;  C)  the  MTF  of  the  system  with  (1)  32°  angular  range,  17  projections 
(blue)  and  9  projections  (red),  and  (2)  32°  angular  range,  17  projections  (blue)  and  9  projections  (red);  D)  Reduction 
in  MTF  as  a  function  of  the  central  scanning  angle;  E)  Confirmation  of  MTF  results  using  a  clinical  tomosynthesis 
system;  F)  Simulated  (red)  vs  actual  (black)  MTF  data. 
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The  effect  of  the  number  of  projections  used  to  create  a  reconstruction  on  the  MTF  was  assessed 
by  reducing  the  number  of  projections  to  a  subset  of  the  17  acquired  projections.  These  data  are 
shown  in  Figure  7A,  and  1C.  When  the  number  of  projections  is  changed  (but  the  total  angle 
spanned  is  held  constant)  there  is  little  effect  on  spatial  resolution.  As  the  number  of  projections 
and  total  angle  spanned  is  changed  together,  there  is  a  significant  reduction  improvement  in 
resolution. 

Separate  validation  of  these  data  was  perfonned  on  a  clinical  tomosynthesis  system  (Figure  7E), 
and  a  simulation  of  the  spatial  resolution  was  generated.  This  simulation  shows  good  agreement 
with  experimental  results. 

The  paper  on  spatial  resolution  (See  Appendix)  is  currently  under  review.  The  paper  on  the 
design  of  the  system  and  the  temporal  resolution  is  currently  being  written. 


4)  KEY  RESEARCH  ACCOMPLISHMENTS: 

Designed  prototype  dynamic  contrast  enhanced  imaging  system  intended  for  proof-of- 
principle  experiments. 

Fabricated  and  assembled  the  components,  including  main  gantry,  filter  wheel, 
collimator,  phantom  holder  and  detector  assembly. 

Fabricated  x-ray  generator  interface,  including  dual-energy  control 

Wrote  all  image  acquisition  software,  including  integration  of  the  x-ray  generator,  motor 

controllers  and  x-ray  detector  software 

Developed  a  method  for  quantitative  reconstruction  of  contrast-enhanced  images  with 
scatter  and  thickness  suppression. 

Designed  a  number  of  key  phantoms 

Performed  experiments  in  quantifying  the  contrast  agent  uptake. 

Performed  experiments  detennining  the  temporal  characteristics  of  the  system 
Performed  experiments  determining  the  spatial  characteristics  of  the  system 
Demonstrated  that  4D  time-resolved  contrast-enhanced  DBT  is  feasible 
Working  to  report  these  results 
Working  to  write  next  grant  to  continue  this  research. 


5)  CONCLUSIONS: 

We  have  successfully  demonstrated  in  proof-of-principle  that  4D  DCE-DBT  can  be  achieved  in  a 
dose-neutral  fashion.  The  temporal  and  spatial  resolution  properties  are  well  understood.  These 
results  will  allow  us  to  write  and  submit  the  next  phase  of  this  research  -  the  construction  of  a 
clinical  system  for  clinical  trial.  In  addition,  we  have  begun  to  write  up  these  results  for 
publication  in  the  peer-reviewed  literature,  and  we  will  seek  patent  protection  for  the  acquisition, 
reconstruction  and  image  processing  methods. 
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ABSTRACT 


PURPOSE 

35  The  purpose  of  this  work  is  to  characterize  the  spatial  dependence  of  the  in-plane 
modulation  transfer  function  (MTF)  in  tomosynthesis  and  investigate  the  influence  of  the 
acquisition  geometry. 

METHODS 

40  The  in-plane  MTF  was  calculated  from  measured  edge  spread  functions  using  a  custom 
test  tool  consisting  of  a  0.250  mm  thick  lead  sheet  affixed  to  a  5.88  mm  acrylic  sheet. 
Images  were  obtained  using  a  non-clinical  bench-top  tomosynthesis  system; 
reconstructions  were  performed  using  customized  commercial  software.  The  dependence 
of  the  MTF  was  measured  with  respect  to  the  following  factors:  1)  the  offset  of  the 
45  reconstruction  plane  from  the  plane  containing  the  edge;  2)  the  obliquity  of  the 
reconstruction  plane  with  respect  to  the  detector;  3)  the  number  of  projections  per 
reconstruction;  and  4)  the  acquisition's  angular  range.  Findings  were  validated  by 
developing  an  analytical  model  of  the  MTF  (accounting  for  blurring  due  to  focal  spot 
size,  magnification,  reconstruction  pixel  size,  out-of-focus  plane,  and  reconstruction 
50  blurring)  and  by  replicating  the  experiment  on  a  clinical  tomosynthesis  system. 

RESULTS 

The  MTF  is  degraded  when  the  edge  is  located  between  two  reconstruction  planes.  The 
degree  of  degradation  of  the  MTF  increases  as  the  distance  between  the  edge  and  the 
55  reconstruction  plane  increases  and  as  the  angular  range  of  acquisition  increases. 
Degradation  increases  very  slightly  as  the  obliquity  of  the  reconstruction  plane  with 
respect  to  the  detector  increases  over  the  range  0°  to  32°.  The  number  of  projections  per 
reconstruction  does  not  affect  the  degradation.  The  simulated  MTFs  generated  by  the 
analytical  model  were  concordant  with  findings  of  the  MTF  degradation  factors,  and  a 
60  validation  experiment  on  a  clinical  tomosynthesis  system  provided  independent 
confirmation  that  the  MTF  degrades  with  increasing  distance  between  the  edge  and  the 
reconstruction  plane. 

CONCLUSION 

65  The  MTF  of  tomosynthesis  systems  is  spatially  anisotropic  and  varies  with  sub-slice 
spacing.  Reconstruction  of  oblique  planes  results  in  minimal  degradation  of  the  MTF  and 
thus  may  be  clinically  acceptable. 


70  Key  words:  Digital  breast  tomosynthesis  (DBT),  4D  imaging,  modulation  transfer 
function  (MTF),  spatial  resolution. 


I.  INTRODUCTION 


Digital  breast  tomosynthesis  has  become  widely  adopted  in  breast  cancer  screening  by 
alleviating  the  obscuration  of  overlapping  breast  features.  Although  tomosynthesis 
75  provides  accurate  three-dimensional  localization  and  visualization  of  features,  its 
effectiveness  is  still  limited  in  the  detection  of  tumors  that  do  not  demonstrate  readily 
discernible  mammographic  features1.  Breast  magnetic  resonance  (MR)  imaging  is 
sensitive  to  such  tumors  because  abnormal  metabolic  activity  associated  with  tumor 
growth  can  be  visualized  using  vascular  contrast  agents,  but  MR  lacks  the  ability  to 
80  visualize  microcalcifications  and  other  subtle  signs  of  many  breast  cancers'”  . 

By  combining  the  merits  of  these  methods,  four-dimensional  dynamic  contrast-enhanced 
(4D-DCE)  tomosynthesis  may  be  able  to  improve  the  accuracy  and  precision  of  breast 
cancer  screening.  4D-DCE  tomosynthesis  involves  the  injection  of  a  radio-opaque 
contrast  agent  through  the  tissue  of  interest,  and  acquiring  sets  of  projections  (each  of 
85  which  produces  a  3D  reconstruction)  at  different  time-points  as  the  contrast  agent  moves 
through  the  tissue.  The  movement  of  the  contrast  agent  through  the  tissue  should  be 
observable  from  analysis  of  the  3D  volume  created  at  each  time-point.  Thus,  4D-DCE 
imaging  would  allow  the  3D  visualization  of  blood  flow  in  the  breast  over  time;  perfusion 
data  can  help  draw  conclusions  about  abnormal  angiogenesis,  which  is  a  well- 
90  characterized  precursor  to  tumor  development5.  At  the  same  time,  the  tomosynthesis  data 
can  provide  accurate  anatomical  information. 
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Obtaining  time-resolved  reconstructions  by  current  tomosynthesis  techniques  would 
require  the  acquisition  of  multiple  series  of  tomosynthesis  data  over  a  number  of  time- 
points,  with  each  reconstruction  produced  by  a  full  set  of  projections  from  each  time- 
95  point.  This  method  would  greatly  increase  the  radiation  dose  to  the  patient.  To  address 
this  issue,  we  have  investigated  a  method  to  generate  four-dimensional  reconstructions 
from  a  limited  set  of  projections  in  which  only  one  projection  is  acquired  per  time -point 
and  each  sequential  projection  is  acquired  at  a  different  angle.  We  hypothesize  that 
accurate  time-series  data  could  be  produced  by  reconstructing  limited  subsets  of  these 
100  projections  (e.g.  applying  a  sliding  window  over  the  set  of  projections),  and  even  using 
subsets  that  do  not  span  the  entire  angular  range  of  acquisition.  Figure  1  shows  a 
schematic  representation  of  a  hypothetical  reconstruction  scheme  in  which  projections  at 
eight  angles  are  acquired  at  different  times,  and  subsets  of  three  projections  are  used  for 
reconstruction. 


Recon  1  Recon  3  .  Recon  7 

Recon  2  Recon  4  Recon  8 


105 

Figure  1:  A  hypothetical  reconstruction  scheme  in  which  projections  at  eight  different 
angles  are  acquired  at  different  time -points  and  two  sweeps  are  made  across  these  eight 
angles.  Three  projections  are  used  per  reconstruction. 
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110  In  order  to  identify  the  lower  limit  of  resolvable  features  of  reconstructions  produced  by 
such  a  4D-DCE  tomosynthesis  system,  we  have  sought  to  characterize  the  spatial 
dependence  of  the  in-plane  modulation  transfer  function  (MTF)  in  tomosynthesis  and 
investigate  the  influence  of  the  acquisition  geometry,  particularly  in  reconstructions  of 
planes  oblique  to  the  plane  of  the  detector.  The  study  of  the  spatial  resolution  of  oblique 
115  reconstructions  is  relevant  to  this  application  because  the  reconstruction  of  3D  images 
from  a  subset  of  projections  may  utilize  a  set  of  projections  acquired  at  angles  that  are  not 
symmetric  about  the  source-detector  axis.  Additionally,  oblique  reconstructions  may 
generally  provide  separate  advantages  for  improving  the  visualization  of  features,  such  as 
microcalcifications. 

120  A  custom  bench-top  tomosynthesis  system  was  used  to  acquire  projections  of  a  custom 
lead  edge  test-tool,  which  were  then  used  to  produce  3D  reconstructions.  The  in-plane 
MTF  was  calculated  from  measured  edge  spread  functions  and  the  dependence  of  the 
MTF  was  measured  with  respect  to  the  following  factors:  1)  the  offset  of  the 
reconstruction  plane  from  the  plane  containing  the  edge;  2)  the  obliquity  of  the 
125  reconstruction  plane  with  respect  to  the  detector;  3)  the  number  of  projections  per 
reconstruction;  and  4)  the  acquisition's  angular  range.  Findings  were  validated  by 
developing  an  analytical  model  of  the  MTF  (accounting  for  blurring  due  to  focal  spot 
size,  magnification,  reconstruction  pixel  size,  out-of-focus  plane,  and  reconstruction 
blurring)  and  by  replicating  the  experiment  on  a  clinical  tomosynthesis  system. 
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130  II.  METHODS 


A.  Image  acquisition 

Projection  images  were  acquired  using  a  custom  x-ray  imaging  testing  station  with  x-ray 
source,  detector,  rotary  stage,  and  filter  as  shown  in  Figure  2.  The  testing  station  uses  a 
GE  Senographe  500  x-ray  generator  (General  Electric  Health  Care,  Chalfont-St.  Giles, 
135  UK)  with  a  selenium  flat-panel  detector  (Direct  Radiography  Corp,  Newark,  DE).  Parker 
Motion  stepper  motors  (Parker  Hannifin,  Rohnert  Park,  CA)  are  used  to  drive  a  rotary 
stage  on  which  a  phantom  can  be  mounted  and  rotated,  as  well  as  a  custom  filter  wheel; 
magnetic  sensors  are  used  to  precisely  control  the  position  of  the  rotary  stage  and  filter 
wheel.  The  filter  wheel  is  laser-cut  acrylic  with  various  mounted  filter  metals,  including 
140  0.5  mm  of  aluminum  used  in  this  experiment.  Additionally,  the  x-ray  tube  has  a  25  pm 

molybdenum  filter  mounted  to  the  exit  port. 

For  the  purpose  of  describing  our  methods,  we  define  an  origin  axis  as  the  line  between 
the  center  of  the  focal  spot  and  the  detector  that  intersects  the  detector  at  a  right  angle.  In 
this  testing  station,  the  source  and  detector  are  stationary  and  the  test  object  moves  by 
145  rotation  of  the  rotary  stage  (the  rotary  stage  is  centered  on  the  origin  axis  such  that 
rotation  of  the  test  object  is  equivalent  to  rotation  of  the  x-ray  tube  and  detector  around 
the  axis  of  the  test  object).  This  simulates  an  acquisition  procedure  in  which  the  source 
and  the  detector  move  around  the  subject;  for  the  purposes  of  describing  our  methods,  the 
term  "angular  range"  will  be  used  to  describe  the  rotation  of  the  test  object  during  image 
150  acquisition.  The  components  of  the  testing  station  are  coordinated  by  software  written  in 
C  (Visual  Studio  2010,  Microsoft,  Redmond,  WA)  which  allows  the  user  to  control  the 
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x-ray  exposure  parameters,  filter  wheel  setting,  acquisition  geometry  parameters,  detector 


acquisition  time,  and  the  acquisition  itself. 


38.5  cm  87.5  cm 


Flat-panel  test  tool 

detector 


Detector  Rotary  stage 

and  test  tool 


Collimator  and  X-ray  tube 
filter  wheel 


155  Figure  2:  (top)  The  upper  image  depicts  a  schematic  of  custom  x-ray  imaging  station 
from  aerial  view.  An  x-ray  source,  aluminum  filter,  and  rotary  stage  with  test  object  are 
positioned  in  succession  in  front  of  a  flat  panel  detector,  (bottom)  The  image  below  is  a 
photograph  of  the  x-ray  testing  station  (filter  wheel  is  not  visible). 

Projections  were  acquired  using  32  kV  and  83  mAs  exposures  through  an  aluminum 
160  filter  of  0.5  mm  thickness  and  a  molybdenum  filter  of  25  pm.  Unless  otherwise  stated,  17 
projections  were  acquired  with  an  angular  range  of  ±16°,  each  at  an  interval  of  2°.  The 
exposure  was  collimated  to  the  size  of  the  object  being  imaged,  in  order  to  reduce  scatter. 
Backprojection  filtering  was  used  to  reconstruct  the  images  (Briona,  Real  Time 
Tomography  LLC,  Villanova,  PA). 
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Figure  3:  (left  top)  Schematic  of  lead  edge  test  tool  used  for  assessment  of  spatial 
resolution.  The  grey  region  is  a  sample  region  of  interest  used  for  calculation  of  the  edge 
spread  function,  (right)  Sample  image  of  reconstruction  in  the  plane  of  the  lead  edge,  (left 

bottom)  Photo  of  lead  edge  test  tool. 

170  A  test  tool  consisting  of  a  0.250  mm  thick  lead  edge  affixed  to  a  5.88  mm  acrylic  sheet 
was  imaged  on  the  system,  designed  as  shown  in  Figure  3.  The  edge  is  tilted  at  a  3°  angle 
and  has  six  small  holes  of  approximately  1  mm  diameter  to  assist  in  obtaining  the  correct 
oblique  reconstruction  plane  for  modulation  transfer  function  analysis  and  to  calibrate  the 
relationship  between  pixels  and  physical  distance  in  the  reconstruction  plane. 

175  B.  Modulation  Transfer  Function  (MTF) 

The  MTF  is  computed  using  the  algorithm  described  by  Carton  et  al.6.  For  the 
tomosynthesis  reconstructions,  the  MTF  is  computed  by  taking  the  MTF  of  an  edge  in  a 
slice  of  the  reconstruction.  Briefly,  the  procedure  is  as  follows: 

1.  A  rectangular  region  of  interest  (of  750  x  350  pixels  or  80.3  x  37.5  mm)  including 

180  the  edge  is  selected. 

2.  A  derivative  filter  is  used  to  find  the  edge. 
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3.  The  edge  is  interpolated  and  a  super-sampled  edge  spread  function  (ESF)  is 
computed  over  the  ROI  using  a  sliding  window  with  a  sub-region  column  size  of 
20. 

4.  A  monotonic  filter  is  used  to  reduce  the  noise  of  the  ESF  . 

5.  The  line  spread  function  (LSF)  is  computed  by  differentiation  of  the  ESF  with  a 
differentiation  kernel  of  [-1  1]  and  correction  for  finite  element  differentiation. 

6.  The  MTF  is  computed  by  taking  the  Fourier  transform  of  the  ESF. 


C.  Experiments 

190  1.  Variation  of  MTF  with  reconstruction  slice  offset 

The  effect  of  the  reconstruction  slice  offset  relative  to  the  plane  that  contains  the  edge 
was  measured  by  obtaining  reconstruction  slices  in  0.1  mm  increments  from  the  true  slice 
depth  containing  the  edge.  The  MTF  of  the  same  region  of  interest  encompassing  part  of 
the  test  tool  as  described  previously  in  section  B  of  the  Methods  was  measured  in  each 
195  image.  The  projections  for  these  reconstructions  are  centered  around  the  origin  axis,  in 
which  the  lead  sheet  is  parallel  to  the  detector  in  the  central  projection  angle  (the  edge 
was  rotated  from  -16°  to  16°  in  2°  increments). 
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Figure  4:  (left)  Graphical  depiction  of  the  process  by  which  the  MTF  is  analyzed  at  each 
200  increasing  reconstruction  slice  offset,  in  which  a  feature  of  interest  is  in  the  smaller  (red) 
plane  and  the  reconstruction  plane  is  shown  as  the  larger  (green)  plane,  (right)  Graphical 
depiction  of  the  intersection  of  the  projections  at  the  plane  containing  a  feature  of  interest 

and  the  reconstruction  plane. 

205 

This  procedure  is  depicted  graphically  in  Figure  4,  in  which  the  gray  cube  represents  a 
volume  of  tomosynthesis  data,  a  feature  of  interest  is  located  in  the  smaller  (red)  plane, 
and  the  larger  (green)  reconstruction  plane  is  offset  from  the  plane  containing  the  feature 
of  interest  by  various  distances  and  the  MTF  measured  at  each  offset.  All  MTF  data  are 
210  corrected  as  will  be  described  in  section  C4  of  the  Methods.  The  diagrams  on  the  right  of 
Figure  4  show  where  the  projections  intersect  the  plane  containing  the  edge  and  the 
reconstruction  plane  (this  diagram  will  be  elaborated  upon  in  section  C4  of  the  Methods). 
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2. 


Variation  of  MTF  with  oblique  reconstruction  angle 


The  effect  of  the  reconstruction  slice  angle  in  an  oblique  reconstruction  on  the  MTF  was 
215  assessed  by  changing  the  angle  of  the  lead  sheet  with  respect  to  the  detector  in  the 
reconstruction.  This  was  done  by  obtaining  projections  of  the  test  tool  in  which  the 
central  projection  was  not  one  in  which  the  test  tool  was  parallel  to  the  detector.  The 
MTF  was  computed  for  the  same  ROI  at  different  depths  for  each  obliquely  angled 
reconstruction.  The  MTF  for  each  reconstruction  slice  angle  was  taken  as  the  curve  with 
220  the  greatest  overall  MTF  for  each  of  the  different  depths  at  that  angle.  This  maximum 
MTF  of  a  group  of  MTFs  was  computed  by  smoothing  the  data  and  looking  for  the  curve 
with  the  highest  frequency  at  which  the  MTF  drops  to  50%.  All  MTF  data  are  corrected 
as  will  be  described  in  section  C4  of  the  Methods. 

3.  Variation  of  MTF  with  the  number  of  projections 

225  The  effect  of  the  number  of  projections  used  to  create  a  reconstruction  on  the  MTF  was 
assessed  by  reducing  the  number  of  projections  to  a  subset  of  the  17  acquired  projections. 
In  one  experiment,  the  angular  range  of  acquisition  was  maintained  at  ±16°  while  the 
angular  spacing  between  projections  was  increased.  Projection  subsets  of  size  15,  13,  11, 
and  9  were  examined.  In  another  experiment,  the  angular  range  was  reduced  to  ±14°, 
230  ±12°,  ±10°,  and  ±8°  while  maintaining  the  angular  spacing  between  projections  (and  the 

projection  subset  size  reduced  accordingly  for  each  experiment).  All  MTFs  are  corrected 
as  will  be  described  in  section  C4  of  the  Methods. 
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4. 


MTF  Simulation 


We  hypothesized  that  the  decay  of  the  MTF  with  respect  to  reconstruction  slice  offset 
235  results  from  a  combination  of  image  blurring  due  to  the  size  of  the  reconstruction  pixels, 
blurring  due  to  the  offset  of  the  reconstruction  slice  to  the  true  slice  that  contains  the 
edge,  the  blurring  due  to  the  reconstruction  process,  and  the  blurring  due  to  the  focal  spot 
and  geometric  properties  of  the  imaging  system.  It  should  be  noted  that  the  MTF  is 
known  to  exhibit  degradation  due  to  the  oblique  incidence  angle  of  x-rays  on  the 

o  i  c 

240  detector  '  .  However,  because  the  detector  used  in  this  experiment  is  selenium  based,  the 
MTF  degradation  due  to  this  effect  is  assumed  to  be  close  to  zero  and  is  not  modeled16. 
The  combination  of  all  of  the  above  mentioned  blurring  factors  can  be  modeled  as  a 
convolution  of  rect  functions  in  the  spatial  domain,  or  a  multiplication  of  sine  functions 

IV 

in  the  frequency  domain  . 

245  In  the  following  equations,  x  =  position  (mm),/=  spatial  frequency  (mm'1),  p  =  pixel  size 
in  the  reconstruction  (pm),  cl  =  reconstruction  slice  offset  (mm),  9  =  half  of  the  angular 
range  (degrees),  m  =  geometric  magnification,  a  =  width  of  the  Gaussian  focal  spot 
distribution  (mm),  Mp=  MTF  of  a  single  projection  image  parallel  to  the  detector,  Mr  = 
MTF  of  a  reconstruction  slice  where  the  central  angle  has  the  edge  parallel  to  the 
250  detector,  and  is  the  scalar  multiplication  operator  (see  Figure  5  for  a  graphical 
depiction  of  geometric  terms).  The  blurring  due  to  the  reconstruction  slice  offset  can  be 
modeled  as  a  rect  function  with  width  equal  to  the  width  of  the  shadow  cast  by  the  feature 
of  interest  onto  the  reconstruction  plane,  or  the  value  of  2b  in  Figure  5. 
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LSF0, fs„<X)  =  red  (1) 


255  The  blurring  due  to  the  size  of  the  reconstruction  pixels  is  also  modeled  as  a  rect  function 

18  20 

with  width  equal  to  the  pixel  size  in  the  reconstruction  plane  '  . 

LSFdetector(x)  vect  ^  (2) 


The  blurring  due  to  the  reconstruction  process  was  modeled  by  dividing  the  MTF  of  a 
reconstruction  in  the  plane  of  the  edge  at  0  mm  slice  offset  by  the  MTF  of  the  central 
260  projection  image  acquired  in  the  direction  of  the  origin  axis.  The  same  projection  is  used 
for  all  simulated  blurring  due  to  reconstruction. 


MTF 


(/)=  £ 


reconstruction 


(3) 


The  blurring  due  to  the  focal  spot  size  was  modeled  as  described  by  Siewerdsen  and 
Jaffray  in  200021. 

771-1 

265  MTF  focalspot(/)=  (4) 

Taking  the  Fourier  transforms  of  LSF0ffset  and  LSFdetector  (ignoring  the  resulting 
coefficient  in  order  to  normalize  the  maximum  value  to  1)  and  multiplying  all  four  factors 
together  results  in  the  following  equation  for  the  MTF  of  a  specific  reconstruction  slice 
offset  from  the  true  edge  plane  and  a  specific  angular  range: 

270  MTF (/)  =  \sinc(nfp)\  *  \sinc(2nfd  tan(0))|  *  —  *  (e  n*  ™  *a*T)2  (5) 

Mp 
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As  illustrated  in  Figure  5,  we  expect  that  the  blurring  is  due  in  part  to  the  offset  of  the 
reconstruction  plane  from  the  true  plane  containing  the  edge,  which  causes  the 
projections  passing  through  the  edge  to  cast  a  magnified  shadow  onto  the  reconstruction 
plane. 


275 
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Figure  5:  Depiction  of  the  blurring  due  to  an  offset  of  the  reconstruction  plane  from  the 
true  plane  that  contains  the  edge  based  on  the  angular  range  of  acquisition. 

In  the  acquisitions  described  in  sections  Cl,  C2,  and  C3  of  the  Methods,  a  low  frequency 

280  artifact  was  observed  in  the  in-plane  MTF  measured  from  their  reconstructions.  It  is 

hypothesized  that  this  artifact  was  due  to  non-linearity  in  the  detector,  which  is  an  older, 

selenium-based  detector  (the  possible  sources  of  this  artifact  will  be  described  in  more 

detail  in  the  Discussion).  This  artifact  is  depicted  in  Figure  6  below,  where  the  MTF  data 

at  83  mAs  and  0  mm  reconstruction  slice  offset  is  superimposed  against  a  simulated  MTF 

285  curve  at  0  mm  reconstruction  slice  offset  (using  the  equation  described  above)  in  which 

the  artifact  is  absent. 
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Figure  6:  Simulated  MTF  data  at  0  mm  slice  offset  superimposed  with  actual  83  mAs 
MTF  data  at  0  mm  slice  offset.  A  low  frequency  artifact  is  visible  in  the  actual  data,  as 
290  well  as  some  offset  in  the  corresponding  MTF  values. 

The  quotient  of  the  simulated  MTF  with  no  slice  offset  and  the  actual  MTF  with  no  slice 

offset  was  computed,  smoothed  (for  all  values  at  greater  than  6  mm  1  spatial  frequency, 

the  data  was  set  to  1  due  to  the  low  signal  of  the  MTF  data),  and  used  as  a  correction 

295  factor  (see  Figure  7);  assuming  the  correctness  of  the  simulation,  this  quotient  would 

represent  only  the  artifact  and  it  should  be  valid  to  apply  this  correction  factor  to  MTF 

data  at  all  slice  offsets. 
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300  Figure  7:  The  quotient  of  the  simulated  MTF  with  no  slice  offset  and  the  actual  MTF 

with  no  slice  offset. 

The  experimental  MTF  at  all  slice  offsets  was  adjusted  by  multiplication  by  this 
correction  factor  obtained  at  0  mm  offset.  If  the  factors  included  in  the  simulation  are 
305  accurate,  it  can  be  expected  that  this  operation  would  effectively  eliminate  the  low 

frequency  artifact  at  all  slice  offsets,  as  well  as  in  reconstructions  with  different  angular 
ranges. 

5.  Experimental  validation  with  a  clinical  system 

310  In  order  to  validate  the  MTF  measurements  seen  in  reconstructions  at  different  slice 
depths  on  the  RTT  system,  a  similar  experiment  was  conducted  to  observe  the  change  in 
MTF  with  reconstruction  slice  offset  on  a  clinical  Selenia  Dimensions  (Hologic,  Bedford, 
MA)  tomosynthesis  system.  The  Hologic  system  reconstructs  slices  between  the  detector 
and  the  breast  support  in  1  mm  increments,  where  the  detector  and  breast  support  are 
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315  parallel.  In  order  to  replicate  the  effect  of  changing  the  reconstruction  slice  depth  at 
smaller  increments,  a  fixture  was  designed  to  elevate  the  test  tool  between  the  detector 
and  the  breast  support.  The  experimental  setup  is  shown  in  Figure  8. 

H 

Figure  8:  Experimental  setup  to  measure  MTF  in  planes  offset  from  the  plane  of  the  lead 
320  edge  on  a  clinical  Hologic  tomosynthesis  system. 

The  Hologic  tomosynthesis  system  has  an  angular  range  of  ±7.5°  and  acquires  a 
projection  every  1.07°.  The  images  were  produced  using  32  kV  and  6  mAs  for  each 
projection.  The  test  tool  was  moved  in  100  pm  increments  for  2  mm,  in  order  to  traverse 
the  distance  equivalent  to  three  adjacent  reconstruction  slices  (which  are  spaced  in  1  mm 
325  steps).  The  MTF  was  measured  using  the  technique  detailed  in  section  B  of  the  Methods. 
Measurements  were  made  for  several  reconstruction  slices  as  the  test  tool  was  translated 
along  the  origin  axis.  We  expected  to  see  results  concordant  with  the  experiment  detailed 
in  section  C.l  of  the  Methods,  in  which  the  reconstruction  slice  offset  was  increased  in 
100  pm  increments. 

330  III.  RESULTS 

A.  Variation  of  MTF  with  reconstruction  slice  offset 

Figure  9  shows  the  change  in  the  MTF  as  the  reconstruction  slice  offset  is  changed.  Each 
MTF  is  computed  from  a  region  of  interest  of  size  750  pixels  by  350  pixels  (80.3  mm  x 
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37.5  mm).  It  can  be  seen  that  the  MTF  with  the  least  decay  is  from  the  plane  that  contains 
335  the  true  position  of  the  edge  (0  mm  offset).  As  the  reconstruction  plane  deviates  from  the 
true  plane,  the  MTF  decays  increasingly.  The  MTF  decay  appears  to  be  symmetric  with 
respect  to  positive  and  negative  position  offsets  of  the  reconstruction  plane. 


Spatial  Frequency  (mm’1) 


Figure  9:  MTF  versus  frequency  of  reconstruction  slices  of  the  test  tool  in  the  plane  of 
340  the  edge,  at  several  selected  offsets  from  the  plane  containing  the  edge.  As  the  curves 
transition  from  superior  to  inferior,  the  reconstruction  slice  offset  is  increased  from  0  to 
0.5  mm.  It  can  be  observed  that  the  MTF  decays  as  the  reconstruction  plane  is  further 
removed  from  the  plane  containing  the  edge  and  that  offset  in  either  direction  is 

symmetric  in  the  MTF. 

345 


The  spatial  frequency  at  which  the  MTF  at  each  reconstruction  slice  offset  (over  a  range 
of  0-2  mm)  drops  to  50%  of  its  maximum  value  is  plotted  in  Figure  10,  in  order  to  better 
depict  the  rate  of  degradation. 
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350  Figure  10:  The  spatial  frequency  at  which  the  MTF  drops  to  50%  of  its  maximum  value 
versus  the  reconstruction  slice  offset  over  a  0-2  mm  range  of  slice  offsets,  in  0. 1  mm 
increments.  The  pictograph  in  the  upper  right  comer  shows  the  original  MTF  data  for  this 

set  of  reconstruction  slice  offsets. 


355  B.  Variation  of  MTF  with  oblique  reconstruction  angle 

It  was  observed  that  the  MTF  changes  slightly  as  the  obliquity  of  the  reconstruction  plane 
is  increased.  Figure  11  shows  the  distribution  of  MTFs  across  a  64°  range  of 
reconstruction  slice  angles  (-32°  to  32°  with  respect  to  the  plane  parallel  with  the  detector, 
in  2°  increments).  There  exists  a  distinct  trend  in  the  data  showing  that  the  MTF  degrades 
360  as  the  obliquity  of  the  reconstruction  plane  increases,  shown  in  Figure  11  by  the 
transition  from  superior  to  inferior  curves  (blue  to  red).  One  outlier  was  removed  from 
this  dataset. 
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Figure  11:  MTF  versus  frequency  for  a  range  of  reconstruction  slice  angles  in  which  the 
365  lead  edge  is  positioned  in  the  plane  for  each  reconstruction.  The  MTF  curves  transition 
from  blue  to  red  as  the  obliquity  of  the  reconstruction  slice  angle  increases.  A  clear  trend 
is  visible  that  shows  the  MTF  degrading  slightly  as  the  obliquity  increases  (from  the 

superior  to  the  inferior  curves). 

370  To  illustrate  this  trend  further,  the  frequency  at  which  the  MTF  drops  below  50%  of  its 
maximum  value  is  plotted  for  each  of  the  33  different  reconstruction  slice  angles  in 
Figure  12.  The  trend  is  clearly  visible  in  this  figure. 

A  second  order  polynomial  fit  was  plotted  through  this  data  as  shown  in  Figure  12.  In  an 
ideal  system,  the  peak  of  the  hyperbola  would  be  located  at  0°,  but  some  outliers  and 
375  noise  exist  in  the  data.  Figure  12  shows  that  the  MTF  drops  by  4.1%  when  comparing  the 
peak  MTF  to  the  MTF  at  a  32°  obliquity  from  the  peak,  based  on  the  trendline  data. 
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Figure  12:  Polynomial  fit  through  50%  of  maximum  MTF  data  for  each  central  angle. 
The  fitted  data  has  equation  displayed  above  with  an  R“  value  of  0.7818. 

380 

C.  Variation  of  MTF  with  the  number  of  projections 

Based  on  the  analytical  model  described  in  section  C4  of  the  Methods,  we  hypothesized 
that  changing  the  number  of  projections  and  preserving  the  angular  range  of  acquisition 
would  not  have  an  effect  on  the  MTF.  Figure  13  shows  the  MTF  in  the  plane  of  the  test 
385  tool  for  a  varying  number  of  projections  used  to  create  the  reconstruction.  There  is  little 
change  in  the  MTF  as  the  number  of  projections  changes,  and  there  is  no  distinguishable 
trend  so  long  as  the  angular  range  remains  the  same.  This  appears  to  be  true  regardless  of 
which  n  projections  are  used  to  reconstruct  a  subset  of  size  n;  that  is,  given  a  set  of  17 
projections,  a  reconstruction  using  any  5  of  them  will  produce  the  same  MTF  as  if  all  17 
390  projections  were  used  in  the  reconstruction,  provided  the  outer-most  projection  angles  are 
included  in  the  five. 
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Figure  13:  MTF  versus  frequency  for  reconstructions  in  the  plane  of  the  lead  edge 
created  using  varying  numbers  of  projections. 

395 

Based  on  the  analytical  model  described  in  section  4  of  the  Methods,  we  hypothesized 
that  the  MTF  would  degrade  at  a  slower  rate  with  respect  to  the  reconstruction  slice  offset 
if  the  angular  range  was  reduced,  due  to  a  reduced  contribution  to  blurring.  As  the 
400  angular  range  of  acquisition  decreases,  the  MTF  decays  more  slowly  with  respect  to  the 
reconstruction  slice  depth.  Figure  14  shows  the  MTF  as  the  reconstruction  slice  depth 
deviates  from  the  plane  containing  the  edge,  in  a  reconstruction  with  9  projections  and  an 
angular  range  of  ±8°  (50%  of  full  angular  range).  It  can  be  observed  that  the  amplitude  of 
the  MTF  decays  more  slowly  than  in  reconstructions  with  a  larger  angular  range. 
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Figure  14:  MTF  versus  frequency  for  a  reconstruction  with  9  projections  and  an  angular 
range  of  +8°,  varied  over  several  reconstruction  slice  depths.  The  reconstruction  slice 
offset  is  increased  as  the  curves  transition  from  superior  to  inferior. 

It  then  follows  that  the  MTF  degradation  with  respect  to  reconstruction  slice  offset  should 

410  remain  the  same  if  the  angular  range  is  preserved  and  any  number  of  projections  is  used 

in  the  reconstruction.  As  can  be  seen  in  Figure  15,  when  the  angular  range  is  kept 

constant  while  the  number  of  projections  is  reduced,  the  MTF  decays  similarly  with 

reconstruction  depth  as  when  the  full  set  of  projections  is  used.  Figure  15  shows  the 

variation  of  MTF  with  reconstruction  slice  depth  in  a  reconstruction  from  9  projections 

415  with  an  angular  range  of  ±16°,  and  this  figure  is  analogous  to  Figure  8,  in  which  all  17 

projections  are  used  in  the  reconstruction. 
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Figure  15:  MTF  versus  frequency  at  different  reconstruction  slice  depths  with  a 
reconstruction  created  from  9  projections  and  an  angular  range  of  ±16°.  It  can  be 
420  observed  that  the  decay  of  the  MTF  is  very  similar  to  that  seen  in  a  reconstruction  with 
the  same  angular  range  but  more  projections.  The  reconstruction  slice  offset  is  increased 
as  the  curves  transition  from  superior  to  inferior. 


D.  MTF  Simulation 

425  In  order  to  gain  a  better  understanding  of  the  sources  of  blurring  in  the  MTF,  an 
analytical  model  of  the  MTF  was  developed  which  composed  all  of  the  understood 
components  of  blurring  as  described  in  section  4  of  the  Methods.  Figure  16  shows  the 
results  of  this  analytical  model  for  an  angular  range  of  ±16°  and  reconstruction  slice 
offsets  ranging  from  0  to  1.9  mm.  The  analytical  model  shows  an  MTF  similar  to  actual 

430  MTF  data  under  the  given  parameters,  and  with  similar  degradation  properties  as  the 
actual  MTF  data.  Although  the  simulated  data  is  very  well  aligned  with  the  experimental 
data,  the  alignment  is  not  perfect,  owing  to  the  lack  of  full  understanding  of  the 
mechanisms  of  blurring  in  tomosynthesis.  It  should  be  noted  that  experimental  data  at  all 
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depths  in  Figure  16  were  corrected  by  the  same  correction  factor  shown  in  Figure  6. 
435  Despite  this  correction  for  the  artifact  caused  by  hardware  limitations,  some 
misalignment  between  the  simulated  and  the  experimental  data  exists,  likely  due  to 
sources  of  blurring  that  are  not  yet  well  understood. 


Spatial  Frequency  (mm-1) 


Figure  16:  (top)  The  results  of  the  analytical  model  of  the  MTF  (red)  superimposed 
440  against  an  actual  MTF  dataset  acquired  from  the  test  tool  (black),  both  over  a  range  of 
reconstruction  slice  offsets  from  0  to  1.9  mm  in  0.1  mm  steps  with  ±  16°  angular  range, 
(bottom  left)  The  simulation  results  for  this  range  of  slice  offsets,  plotted  alone,  (bottom 
right)  The  experimental  MTF  for  this  range  of  slice  offsets,  plotted  alone.  The  superior 
curves  are  for  a  slice  offset  of  0  mm  and  the  inferior  curves  are  for  a  slice  offset  of  1.9 
445  mm. 

In  addition  to  simulating  the  effect  of  increasing  the  reconstruction  slice  offset  at  ±16° 
angular  range,  a  simulation  of  this  effect  at  +8°  angular  was  also  performed  and  compared 
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to  an  actual  dataset,  illustrated  in  Figure  17.  This  simulation  reflected  the  rate  of 
degradation  of  the  actual  data  nicely,  but  once  again,  the  data  was  not  perfectly  aligned. 
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Figure  17:  (top)  The  results  of  the  analytical  model  of  the  MTF  (red)  superimposed 
against  an  actual  MTF  dataset  acquired  from  the  lead  edge  test  tool  (black),  both  over  a 
range  of  reconstruction  slice  offsets  from  0  to  1.9  mm  with  ±  8°  angular  range,  (bottom 
455  left)  The  simulation  results  for  this  range  of  slice  offsets,  plotted  alone,  (bottom  right) 
The  experimental  MTF  for  this  range  of  slice  offsets,  plotted  alone.  The  superior  curves 
are  for  a  slice  offset  of  0  mm  and  the  inferior  curves  are  for  a  slice  offset  of  1 .9  mm. 


E.  Experimental  validation  with  a  clinical  system 

460  The  same  trend  of  decaying  MTF  with  increasing  reconstruction  slice  offset  is  visible  in  a 
repetition  of  the  experiment  on  a  clinical  Hologic  system.  Figure  18  shows  the  MTF 
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versus  frequency  for  several  instances  of  a  single  reconstruction  slice  in  which  the  test 


tool  is  moved  away  from  the  slice  containing  the  edge  in  100  pm  increments. 


Spatial  Frequency  (mm-1) 


465  Figure  18:  MTF  versus  frequency  in  Hologic  tomosynthesis  reconstructions  where  the 
lead  edge  is  moved  away  from  the  reconstruction  slice  that  contains  the  edge.  The  MTF 
can  be  seen  to  decay  as  the  slice  moves  further  from  the  slice  containing  the  edge  (as  the 
curves  transition  from  superior,  red,  to  inferior,  blue). 

470 

Figure  18  appears  to  show  a  consistent  degradation  in  the  MTF  as  the  offset  of  the 
reconstruction  plane  from  the  plane  containing  the  test  tool  increases.  In  order  to  better 
visualize  the  characteristics  of  the  degradation  with  increasing  slice  offset,  the  frequency 
at  which  the  MTF  drops  below  50%  of  maximum  is  plotted  for  each  reconstruction  slice 
475  offset  in  Figure  19,  showing  a  largely  linear  degradation  over  the  2  mm  range.  Compared 
to  Figure  9,  it  is  hypothesized  that  this  data  is  still  degrading  in  the  linear  range, 
analogous  to  the  0  mm  -  0.25  mm  slice  offset  range  of  Figure  9.  For  clinical  purposes,  the 
Hologic  system  generates  reconstruction  slices  at  1  mm  increments,  and  thus,  the 
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maximum  offset  of  the  closest  reconstruction  plane  to  a  feature  of  interest  never  exceeds 


480  0.5  mm,  at  which  there  is  a  maximum  MTF  degradation  of  only  9%. 


Figure  19:  Reconstruction  slice  offset  versus  spatial  frequency  at  which  the  MTF  drops 
below  50%  of  maximum.  This  data  is  fit  linearly  (with  an  R2  value  of  0.977)  and  appears 
to  follow  such  a  linear  degradation  as  the  reconstruction  plane  moves  away  from  the 
485  plane  containing  the  feature  of  interest  over  a  2  mm  range. 


IV.  DISCUSSION 
A.  Consequences  of  this  work 

1.  Effect  of  slice  alignment  on  the  conspicuity  of  clinical  features 

490  Blurring  of  features  in  a  tomographic  image  is  a  non-trivial  problem  that  is  not  present  in 
mammography  or  stereomammography.  In  the  case  of  mammography,  the  attenuation  of 
a  feature  of  any  depth  or  thickness  is  projected  directly  onto  the  image  plane,  whereas  a 
tomographic  reconstruction  interrogates  a  volume  of  data  in  which  features  are  located  in 
certain  planes  depending  on  the  thickness  of  the  feature.  Features  of  small  thickness 

495  along  the  detector-source  axis  but  with  high  attenuation  which  are  readily  visible  in 
mammographic  images  may  be  completely  absent  from  tomographic  reconstructions  due 
to  the  additional  blurring  factors  associated  with  tomosynthesis.  This  has  a  serious  impact 
on  the  use  of  breast  tomosynthesis  in  a  clinical  setting. 
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Multiple  clinical  studies  have  reported  that  tomosynthesis  is  limited  in  its  ability  to  detect 
500  breast  microcalcifications.  Calcifications  in  breast  tissue  are  indicative  of  local  cell  death 
or  secretions  due  to  pathological  activity  and  can  be  precursors  to  breast  cancers  such  as 
ductal  carcinoma  ’  .  Breast  microcalcifications  are  often  100-300  pm  in  thickness, 

whereas  tomographic  reconstruction  slice  spacings  are  often  placed  at  a  500-1000  pm 
offset.  Given  a  300  pm  calcification  and  a  1000  pm  slice  spacing,  there  is  a  70%  chance 
505  that  a  random  set  of  slices  with  1000  pm  spacing  will  not  have  a  single  slice  that  passes 
through  the  calcification.  The  calcification  may  still  be  visible  in  planes  offset  from  the 
actual  position  of  the  calcification,  but  the  spatial  resolution  will  be  degraded  due  to 
blurring.  The  results  presented  in  this  work  may  help  determine  tolerable  reconstruction 
offsets,  and  hence,  the  minimum  slice  spacing  which  retains  sufficient  spatial  resolution 
510  to  depict  small  clinical  features  such  as  microcalcifications. 

The  sources  of  blurring  in  tomosynthesis  are  not  fully  understood.  It  can  be  interpreted 
from  the  results  of  this  work  that  blurring  is  at  least  partially  due  to  a  combination  of:  1) 
the  offset  of  the  reconstruction  plane  from  the  plane  containing  the  feature;  2)  the  angular 
range  of  the  acquisition;  3)  the  detector  element  and  source  properties;  and  4)  the 
515  acquisition  geometry.  As  illustrated  in  Figure  4,  when  the  offset  of  the  reconstruction 
plane  from  the  feature  increases,  the  feature  projects  a  larger  and  more  blurred  shadow 
onto  the  reconstruction  plane.  This  effect  can  be  characterized  by  the  degradation  of  the 
MTF  of  an  edge  as  the  edge  moves  away  from  the  reconstruction  plane,  as  can  be  seen  in 
Figure  8.  The  degradation  is  quite  severe  in  a  system  which  uses  a  ±16°  angular  range  of 
520  acquisition  as  demonstrated  in  section  A  of  the  Results.  In  reconstructions  with  an 
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angular  range  of  ±16°,  the  frequency  at  which  the  MTF  drops  to  half  of  its  maximum 
value  occurs  at  about  a  600  pm  offset  of  the  reconstruction  plane.  On  the  other  hand, 
when  the  angular  range  is  reduced  to  ±8°,  the  MTF  dropoff  with  slice  offset  is  much  less 
steep,  as  depicted  in  Figure  14. 

525  The  rate  of  increase  of  the  blurring  with  increasing  distance  between  the  feature  and  the 
reconstruction  plane  is  due  in  part  to  the  angular  range  of  acquisition.  When  the  angular 
range  of  an  acquisition  is  increased,  the  tomographic  ability  (or  the  ability  to  resolve 
characteristics  of  features  in  three  dimensions)  of  the  reconstruction  is  increased-  .  As 
illustrated  in  Figure  3,  for  a  single  reconstruction  plane  with  some  offset,  as  the 
530  acquisition's  angular  range  increases,  the  shadow  projected  by  the  feature  onto  the  plane 
widens,  representing  an  increase  in  blurring.  Likewise,  as  illustrated  in  Figures  14  and  15, 
decreasing  the  angular  range  decreases  the  rate  at  which  blurring  increases  due  to  the 
reconstruction  slice  offset.  These  results  are  in  line  with  the  observations  of  Maidment  et 
al.  (2005),  who  noted  that  using  a  smaller  range  of  projection  angles  results  in  less 
535  blurring  of  out  of  plane  structures,  and  who  determined  a  relationship  between  resolution 
at  the  limiting  frequency  and  angular  range  in  tomosynthesis  images-  .  Thus,  increasing 
the  angular  range  increases  the  need  for  accurate  slice  spacing  as  features  will  become 
out-of-focus  even  with  small  changes  in  reconstruction  plane  depth.  Practically,  the 
selection  of  the  angular  range  represents  a  trade-off  between  accurate  localization  of  an 
540  existing  feature  in  the  detector-source  axis  (e.g.  characterization  of  the  size  and  location 
of  an  existing  calcification)  and  ease  of  screening  for  features  in  a  three-dimensional 
space  (e.g.  screening  tomographic  reconstructions  to  search  for  microcalcifications). 
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Interestingly,  the  number  of  projections  per  reconstruction  did  not  affect  the  MTF  so  long 
as  the  reconstruction  plane  offset  and  the  acquisition's  angular  range  were  held  constant, 
545  as  shown  in  Figure  15.  This  knowledge  could  potentially  be  used  to  reduce  the  radiation 
dosage  during  tomosynthesis  if  prior  knowledge  of  the  structure  being  imaged  exists  and 
the  known  features  do  not  require  high  resolution  between  the  sparsely  spaced  views.  For 
example,  this  technique  is  particularly  applicable  for  needle  tracking  during  biopsy 
procedures.  Based  on  the  position  of  the  needle  and  the  tissue,  projection  views  at  certain 
550  time  points  may  have  obstructions  that  block  the  view  of  the  feature  or  may  not  contain 
useful  geometric  information.  Selectively  reconstructing  images  from  certain  projections 
based  on  this  knowledge  could  improve  temporal  resolution  or  reduce  radiation  dose. 
However,  decreasing  the  number  of  projections  per  reconstruction  can  lead  to  obvious 
artifacts  from  high  constrast  features.  Small  numbers  of  projections  with  high-contrast 
555  features  generate  artifacts  in  reconstructions,  which  are  blurred  and  made  less  obvious  as 
the  number  of  projections  increases.  Increasing  the  number  of  projections  per 
reconstruction  distributes  the  intensity  of  high  contrast  features  over  several  components 
and  causes  artifacts  to  be  less  obvious. 

2.  Anisotropy  of  resolution 

560  It  is  well  known  that  tomosynthesis  produces  images  with  anisotropic  spatial  resolution, 
due  to  the  limited  angular  range  of  acquisition26.  The  present  work  further  illustrates  the 
anisotropic  nature  of  spatial  resolution  in  tomosynthesis  as  objects  move  in  and  out  of 
focus.  The  characterization  of  this  blurring  based  on  the  aforementioned  factors  can  allow 
for  the  estimation  of  the  appropriate  reconstruction  slice  spacing  for  a  particular  clinical 
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565  problem.  In  order  to  resolve  small  clinical  features  with  an  acceptable  spatial  resolution, 
finer  reconstruction  slice  spacings  may  be  necessary.  Finer  slice  spacings  would  require 
storage  of  a  larger  amount  of  data  or  use  of  a  dynamic  reconstruction  scheme  in  which 
slices  can  be  recomputed  from  projection  data  dynamically  based  on  a  selected  slice 
spacing27. 

570  3.  Effect  of  obliquity  on  spatial  resolution  in  tomosynthesis 

Acciavatti  et  al.  (2012)  have  previously  shown  that  digital  breast  tomosynthesis  is 
capable  of  obliquely  angled  reconstructions  .  Oblique  reconstructions  allow  the 
visualization  of  surfaces  that  are  not  parallel  to  the  detector,  and  have  great  potential  in 
applications  such  as  real-time  imaging.  For  example,  real-time  three-dimensional 
575  tracking  of  a  biopsy  needle  or  other  surgical  tools  is  currently  performed  by  careful 
insertion  of  the  tool  in  a  plane  that  is  parallel  to  the  detector.  Oblique  reconstruction 
would  allow  the  alignment  of  the  reconstruction  slice  with  the  tool  in  the  image 
reconstruction  procedure  rather  than  by  physical  alignment  of  the  tool  during  surgery. 
Additionally,  oblique  reconstructions  may  improve  visualization  of  microcalcifications 
580  without  increasing  radiation  dose. 

The  present  work  demonstrates  that  spatial  resolution  degrades  minimally  as  the  angular 
position  of  the  reconstruction  plane  deviates  from  the  angular  position  of  the  plane  of  the 
detector.  The  MTF  can  be  used  to  characterize  the  in-plane  spatial  resolution  of  obliquely 
reconstructed  planes.  Figure  12  shows  that  the  dropoff  of  the  frequency  value  at  which 
585  the  MTF  reaches  50%  of  the  maximum  MTF  degrades  at  an  increasing  rate  as  the 
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obliquity  increases.  The  maximum  obliquity  tested  in  this  experiment  was  a  plane  angled 
at  32°  to  the  plane  of  the  detector.  At  this  degree  of  obliquity,  it  was  observed  that  the 
frequency  value  at  which  the  MTF  reaches  50%  of  the  maximum  fell  by  only  about  10%, 
indicating  that  spatial  resolution  is  relatively  well  preserved  in  oblique  reconstructions. 
590  This  result  indicates  that  reconstruction  of  oblique  planes  may  be  clinically  acceptable. 

B.  Limitations  of  this  work 

Reconstructions  in  this  work  were  performed  using  backprojection  filtering.  The 
reconstruction  algorithm  supports  super-resolution;  that  is,  sub-pixel  resolution^  .  Super- 

in 

resolution  has  also  been  demonstrated  to  exist  in  iterative  reconstructions  .  The  spatial 
595  resolution  characteristics  of  reconstructions  computed  by  backprojection  filtering  and 
iterative  reconstruction  may  be  different,  and  the  present  work  has  not  investigated  the 
performance  of  iterative  reconstruction.  Lu  et  al.  demonstrated  the  effect  of  voxel  size  on 

o  1 

the  in-plane  intensity  profile  of  iterative  reconstructions  with  subpixel  resolution  .  Based 
on  the  concordance  of  backprojection  filtering  and  iterative  reconstruction  super- 
600  resolution  results,  it  is  likely  that  the  factors  affecting  the  MTF  described  in  this  work 
also  apply  to  iterative  reconstructions.  However,  exploring  this  relationship  is  outside  of 
the  scope  of  the  present  work. 

Although  not  presented  here,  the  effect  of  the  reconstruction  slice  offset  on  the  MTF  in 
oblique  reconstructions  was  examined  and  was  found  to  be  concordant  with  the  MTF 
605  degradation  in  non-oblique  reconstructions  (Figure  8  and  section  A  of  this  work).  The 
present  work  has  not  investigated  the  effect  of  the  angular  range  of  acquisition  on  the 
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degradation  of  spatial  resolution  with  reconstruction  slice  offset  in  oblique 


reconstructions.  This  is  a  direction  of  interest  that  warrants  future  study. 


The  custom  tomosynthesis  system  on  which  the  experimental  projections  were  acquired 
610  is  a  non-clinical  "sandbox"  system.  It  employs  an  older  selenium  flat-panel  detector 
which  is  limited  in  its  linear  range.  The  MTF  measurements  of  reconstructions  from  this 
system  show  a  sharp  drop  in  the  0-1  mm1,  which  is  an  artifact  of  the  detector  being 
driven  into  nonlinear  operation.  Reducing  the  mAs  during  image  acquisition  can  reduce 
this  artifact  but  not  entirely  remove  it.  Figure  20  shows  the  MTF  with  variation  in 
615  reconstruction  slice  offset  of  0  mm  to  ±1  mm  with  a  6  mAs  acquisition  (analogous  to 
Figure  7).  The  intensity  profile  of  the  entire  width  of  the  test  tool  of  the  phantom  was 
symmetric  and  reducing  the  field  size  by  collimation  did  not  affect  the  intensity  profile, 
which  ruled  out  scatter  from  the  acrylic  plate  supporting  the  lead  edge  as  only  one  side  of 
the  lead  edge  was  adjacent  to  the  acrylic. 
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Figure  20:  (left)  Comparison  of  MTF  between  80  mAs  acquisition  and  6  mAs 
acquisition.  The  low  mAs  acquisition  shows  a  much  smaller  non-linearity  artifact  at  low 
frequency,  suggesting  that  this  artifact  may  be  due  to  detector  limitations,  (right) 
Example  of  MTF  data  with  variation  in  reconstruction  slice  offset  at  low  mAs  at 
625  reconstruction  slice  offset  of  0  mm  to  ±1  mm. 
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V.  CONCLUSION 


The  MTF  of  tomosynthesis  systems  is  spatially  anisotropic  and  varies  with  sub-slice 
spacing.  Angular  range  and  slice  alignment  both  contribute  to  the  blurring  that  causes 
630  degradation  of  the  in-plane  MTF.  This  blurring,  due  to  sub-slice  misalignment  of  objects, 
may  adversely  affect  the  conspicuity  of  small  clinical  features  such  as  microcalcifications 
and  may  require  finer  slice  spacing  and  the  storage  of  an  increased  amount  of  data.  An 
increase  in  the  angular  range  of  acquisition  also  contributes  to  blurring  but  improves 
tomographic  ability.  Additionally,  reconstruction  of  oblique  planes  results  in  minimal 
635  degradation  of  the  MTF  and  thus  may  be  clinically  acceptable.  It  is  evident  that  many 
factors  affect  blurring  and  image  quality  and  that  careful  deliberation  on  each  factor  and 
tradeoff  may  be  necessary  to  best  identify  certain  clinical  features  in  both  3D  and  4D 
tomosynthesis. 
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